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Abstract - We study the out of equilibrium steady state properties of the Bose-Fermi-Kondo 
model, describing a local magnetic moment coupled to two ferromagnetic leads that support 
bosonic (magnons) and fermionic (Stoner continuum electrons) low energy excitations. This model 
describes the destruction of the Kondo effect as the coupling to the bosons is increased. Its phase 
diagram comprises three non-trivial fixed points. Using a dynamical large- A/" approach on the 
Keldysh contour, we study two different non-equilibrium setups: (a) a finite bias voltage and (b) 
a finite temperature gradient, imposed across the leads. The scaling behavior of the charge and 
energy currents is identified and characterized for all fixed points. We report the existence of a 
fixed-point-dependent effective temperature, defined though the fluctuation dissipation relations 
of the local spin-susceptibility in the scaling regime, which permits to recover the equilibrium 
behavior of both dynamical and static spin susceptibilities. 
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Introduction. — Understanding the physical proper- 
ties of correlated systems away from thermal equilibrium 
has been a subject of intense research. The current inter- 
est in the theoretical description of far-from-equilibrium 
dynamics is partly driven by recent experimental achieve- 
ments to probe non-thermal states in a controlled fashion 
[l]-[3]. Also, the ability to numerically address the time 
evolution of relatively large systems in conjunction with 
exact methods has lead to new insights, e.g. a better 
understanding of thermalization properties of interacting 
systems At present, it appears that the interplay be- 
tween correlation effects and non-thermal boundary con- 
ditions can lead to a plethora of possible behaviors j^-^. 
The existence of multiple steady-states selectively chosen 
by particular initial conditions [7 and of recurrent non- 
equilibrium solutions ^ are examples of such rich phe- 
nomena. Already classical systems away from equilibrium 
can show rich behaviors. A particularly interesting con- 
cept that emerged from the study of classical scale-free sys- 
tems is the notion of effective temperatures based on local 
non-equilibrium fluctuation-dissipation relations 9pQ . In 
the context of scale-free quantum systems, the concept of 
effective temperatures has so far received only limited at- 
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Quantum critical systems, the quantum- 
analogs of classical scale-free systems, should be sensitive 
to any out-of-equilibrium drive, due to their gapless, scale- 
invariant spectrum. Unlike their classical counterparts, 
quantum critical systems link dynamical and static prop- 
erties already at the equilibrium level. This raises the 
question if the concept of effective temperatures can be 
carried over to quantum scale-free systems in a meaning- 
ful way. While no such well-defined effective temperature 
was found in the ohmic spin-boson model 11 , a prelimi- 



nary study for a quantum critical system concluded that 
the notion of effective temperature based on an extension 
of the fluctuation-dissipation theorem can be meaning- 
ful 
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For quantum critical systems possessing a gravity 
dual, it was recently suggested that an effective temper- 
ature characterizing the out-of-equilibrium current noise 
naturally emerges from the holographic mapping flS^. 

Studying the quantum relaxational {uj < T) regime near 
quantum criticality faces methodological difficulties al- 
ready at the equilibrium level and only few exactly solvable 
cases are available. Extending the study to non-thermal 
boundary conditions, poses further challenges. 

In this letter, we address the out-of-equilibrium prop- 
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erties of current carrying steady states near Kondo- 
destroying quantum criticality in the spin-isotropic Bose- 



Fermi Kondo model (BFKM) 17 -19 consisting of a quan- 
tum spin coupled to a fermionic and a sub-ohmic bosonic 
bath 20,21 . We address this model in terms of a dynam- 
ical large-N limit which gives access to the full quantum 
relaxational dynamics and treats equilibrium and out-of- 
equilibrium correlations on the same footing. This allows 
for a controlled comparison of the fluctuation-dissipation 
theorem with its out-of-equilibrium counterpart. 

Model. — We consider the \aige-N version of the 
multi-channel BFKM (sketched in Fig{l]-(a)) where the 
spin degree of freedom {S) is generalized from SU{2) 
to SU{N) [20)[22]. The fermionic excitations (c) of the 
bath transform under the fundamental representation of 
SU{N) X SU{M) with TV the spin and M the charge chan- 
nels. An out-of-equilibrium steady state is obtained by 
coupling the dot to two sets of such baths {i.e. leads) kept 
at different thermodynamic potentials 6 . The system is 
described by the Hamiltonian: 
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Ho^Hj, (1) 



ql 



(3) 



where a and a are respectively the SU{N)-spm and 
5'/7(M)-channel indices, I = L^R labels the left and 
right leads and p^q are momentum indices. The co- 
tunneling terms in Eq.([2| contain the local operators 

^Iw = ^ Epp'aa' ^aaita'Cp'cr'aU [23 , with t the funda- 
mental representation of SU{N). In terms of their mo- 
mentum counterparts, the local bosonic fields write as 

= ^T^q^i (with i = l,...,Ar). Here, and 
are the number of fermionic and bosonic single particle 
states taken to be proportional to the volume of the leads 
and set to infinity at the end of the calculation. 

The impurity's 5'[/(A^) -operators can be written, in 
terms of pseudo-fermions, S'^ = ^^^ra' fa^a a' f(^' ^ where 
i runs over the N'^ — 1 generators of SU{N) and r is an 
anti-symmetric representation of SU{N) fixed by impos- 
ing the constraint Q = "^^a fafo- = the total number 
pseudo-fermions. In the path integral formalism, this con- 
straint is enforced by a dynamical Lagrange multiplier A. 

Assuming that the exchange interaction derives from an 
Anderson-like impurity model [2^ the matrix Ju^ is given 
by Jw = yfJiJv where Ji ex if/U., with ti the hopping 
term between the level and the leads and U the repulsive 
energy at the impurity site. 

The spectral function of the bosonic bath is 
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Fig. 1: (a) A sketch of the non-equilibrium setup of the BFK 
model, (b) The phase diagram of the BFK model encompass- 
ing three fixed points: multichannel Kondo (MK), critical (C) 
and local moment (LM) . The parameters J and g are the 
couplings of the local moment to the fermions and bosons, re- 
spectively. The arrows denote the RG flow to the three fixed 
points, over-screened Kondo, LM and critical, (c) Sketch of 
the Keldysh contour 7 = 7+ + 7_ . (d) Large- A/" self-energy 
diagrams corresponding to Eqs. (5j6). 
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for \uj\ < A, where A is some cut-off scale. The fermionic 
bath is characterized by a constant density of states at the 

Fermi level pc{oj = 0) that for ^ = yields to the Kondo 

1 

energy scale Tk = Pc{0)~ e (Jl+Jr)pc(o) ^ 

The phase space of this model encompasses three fixed 

points, located in the T = 0, V = hyperplane, that can 

be accessed by varying the ratio g/J^ see Fig{l}(b). An 

over-screened multichannel Kondo phase (MK) at small 

g/J is separated from a local moment (LM) phase by an 

unstable critical point (C). The characterization of each 

phase and the scaling laws for the different quantities were 

obtained in | p!9{[2Q] both by perturbative RG and large- 

methods. 

We consider a non-equilibrium setup where the two 
leads, initially decoupled from the impurity (for t < to), 
are held at chemical potentials Pl = —pn = |e| "1^/2 and 
at temperatures Tl and Tr. As the leads are considered 
to be infinite reservoirs, its bosonic and fermionic distribu- 
tions functions are given respectively by n^^i {00) = 

and Uf^iiuj) = with pi = T-\ At t = to 

the coupling between the leads and the impurity (Hj) is 
turned on. We address the steady state regime by formally 
setting to —00. 

Dynamical large-N out-of-equilibrium. We follow the 
dynamic large- approach of Ref. [20] , generalized to an 
out-of-equilibrium setup in Ref. p]. The generating func- 
tion writes Z = J DcDd'D^D^''DfDpD\e-^, with S 
the action associated with Eq.([2| and where the integra- 
tion of the fields is performed over the forward (7+) and 
backward (7-) Keldysh branches (see Figjl]-(c)). Our ap- 
proach generalizes the calculations under the equilibrium 
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conditions f2Q[|22| to the Keldysh contour Q [24] , 

In the steady-state regime Ga t') = Ga{t — t') (with 
a = f,B) and the Keldysh equation G>'< = GfS>'<G^ 
holds for the local pseudo-particle Green's functions. The 
saddle-point equations simplify to 



it) = 



S>'< it) 



-iG< it) 



iG)'< {t)G<^>{-t) 
-iKG>'<{t)G>^<{t) 
+ig'G>'< (t) [G>'< it) 

q, 



G 



(5) 
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where q = k = = gLQL+gRgR, and J 

The local quantities Gc {t) = 
7^ 7 [^lGcpl (t) + JrGc^pR (t)] and G$ {t) 
^ Eg ^ [ghghG^qL (t) + gRgRG^qR (t)] are defined in 
terms of the lead's fermionic and bosonic Green's functions 
and the Green's functions of the local fields are given by 
GJ^ {t) = -{dt-fi)-^f (t) and G^^ {t) = J'^ - (t), 
with —iji the saddle point value of the Lagrange multiplier 
field A. 

For the numerical evaluation of the saddle point 
equations, we use pc\^) = -We^voy 



OH 



D is Si hight en- 



tD 

where 



Ar(^) 

ergy cutoff of the fermionic density of states, A = 0.1 
hi = 1/2 and q = 1/2. For definiteness we also take 
= 5/4. 

The self-consistent equations were solved iteratively on 
a logarithmic discretized grid with 250 points ranging 
from —lOD to lOD. The criterium for convergence was 
that the relative difference of two consecutive iterations 
was less than 10~^. We also set Jpc (0) = O.S/tt, yielding 
a critical value of the bosonic coupling of gc — 0.2199J. 
We checked that our results hold for other choices of pa- 
rameters. 

The scaling properties of the model at equilibrium are 
analyzed in Ref. 
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The scaling exponents of the 
auxiliary Green functions Gf and Gb are summarized in 
table [l] where a / and a b are defined through 



1-aa 



<Saioo/T) {a = f,B), 



(8) 



and ^(x) is a smooth scaling function with ^(0) ^ 0. 
The imaginary part of the spin susceptibility behaves as 



where 



2a/ - 1. 



(X 



sgn (cj) |cj|' 



(9) 



-•^The SU (M) symmetry of the model ensures that the propagator 
of the local bosonic field B is independent of the channel index a. 
The singular nature of the matrix Jnf substantially simplifies the 
treatment as the combination a/ Jr / JB l — \/ JlI J Br decouples 
from the equations and one is left with a single scalar field B = 
^/Jl/JBl + ^/Jr/JBr. 

= 1/2 corresponds to a particle-hole symmetric setup. Note, 
the out-of-equilibrium conditions considered here respect particle- 
hole symmetry. 
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Table 1: Scaling exponents of the auxiliary Green functions 
Gf and Gb here evaluated for k = 1/2 and = 5/4. 

We report results for two out-of equilibrium situations: 
(a) a finite bias voltage applied across the leads jul = 
—jj^R = V/2 kept at the same temperature T and (b) a 
finite temperature gradient AT = Tl — Tr with fiL = 

llR = 0. 

Observables. The particle number and energy cur- 
rents, denoted respectively J7n and J7e, obtained from the 
continuity equation, are Jb,a,i = —dt {Qb,a,i (^)), where 
b = n^e (for particle and energy current respectively), a = 
c, ^ (for the fermionic and bosonic fields), Qn,a,i = ^a,i is 
the number of particles a on the lead / and Qe,a,i is the 
part of the Hamiltonian of the particles a restricted to the 
lead I. Explicitly, one gets 



Jb. 



^,1 



-^77=^ ^^^^^^ [d^^'cra' (PaU'^qil) " Cc] , 



Jb,c,l 



N 
1 

'n 



OiCTcr' \kp 

ft f .J 
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where V 

Jn,c,p 



7^ 



, Jn,^,q — Jn,c,p — 1? Je,^,q — , 

All operators inside brackets are computed 
at equal time. Note that for the fermions the identity 
sJb,c,L = —sJb,c,R immediately follows. This does not 
hold for the bosons as their number is not conserved by 
the Hamiltonian. In terms of the particle current, the 
total electric current leaving the left lead is given by 
I = —\e\ Jn,c,L, where |e| is the absolute value of the 
electron charge. 

The correlation functions appearing in m- 
volve impurity as well as bath degrees of freedom. As 
Wick's theorem no longer holds, the approach taken here 
consists in inserting external sources, conjugate to the c, <l> 
and / fields, that respect the Keldysh structure. By vary- 
ing the action with respect to the sources one obtains the 
desired correlation functions. The computation of the ex- 
plicit form of generic observables in terms of the distri- 
bution functions of the local and bath fields turns out to 
be rather involved. This cumbersome approach is neces- 
sary as T-matrix-based arguments are not ensured to hold 
for the large-TV generalization of the BFKM. Interestingly 
our results are qualitatively compatible with the T-matrix 
based derivations [6,,20^ whenever a comparison is possible. 

Effective temperature. The steady-state fiuctuation 
dissipation ratio (FDR) is defined, for dynamical observ- 
able O (t,t') = 0(t- t'), as FDRo {u) = [^^(^)+^^('") 



(11) 



and 



0>iuj)-0<iuj) 
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with 0^/< {(jj) being the 
greater/lesser components: 



Fourier-transforms of the 
0>/<{t,t') for t e 7_/+ 



and t' G 7+/- (the definition of the Keldysh branches 
7-1- is given in Figjlj-c). At equilibrium, the fluctua- 
tion dissipation theorem fixes FDRo {oj) = FDReq {oj) = 
tanh {l3uj/2)~^ (with = ±1 for bosonic or fermionic oper- 
ators) uniquely. For a generic out-of-equilibrium system, 
the functional form of the FDR differs from the equilib- 
rium one. 

An observable and frequency dependent "effective tem- 
perature", Pef[,o{^)^ can be defined by requiring that 
tanh(/3eff,o (^) ^/2)"^ FDRq (cj) For the 

regimes where the out-of-equilibrium drive is much smaller 
then the temperature, linear response theory can be used 
and the equilibrium functional form is expected to hold. 
Here we follow Refs. [9j[ll] and define an effective temper- 
ature for the observable O by its asymptotic low frequency 
behavior 

-C 



/^eff,o = /^eff,o (0) lim 



FDRo {(^y 
ujj2 



(12) 



As shown below the effective temperature defined in this 
way holds for all frequencies in the scaling regime. From 
this expression, one obtains an observable-dependent effec- 
tive temperature. Here we will consider, Teffj,Teff,B and 
^eff,x 5 computed from Eq.(12) using the pseudo-particle 



Green's functions and the impurity susceptibility. A pre- 
liminary study of the spin susceptibility near the quantum 
critical LM regime indicated that an effective temperature 



can be defined in the corresponding scaling regime 13 
In this letter, we demonstrate that the notion of effective 
temperature can be successfully used whenever the system 
displays critical scaling in the nonequilibrium regime. The 
evaluation of the pseudo-particle quantities has two pur- 
poses: (i) First, the agreement between these quantities 
indicates that it is possible to define an observable inde- 
pendent effective temperature, (ii) Secondly, as Wick's 
theorem applies to higher spin correlators in the large- 
limit, such observables will automatically have the same 
effective temperature Teff. 

Results . — 

(a) Finite bias voltage. Fig j2] shows the behavior of 
the spin-susceptibility of the impurity and the effective 
temperatures at the MK (^ = 0), C {g = gc) and LM 
{g = 2.5gc) fixed points for F 7^ and AT = 0. 

The effective temperatures (Eq.([T2|), given in Figj2]- 
(right panel), are obtained for T = 5 x IO'^Tk by varymg 
V alone. The full scaling form of Tq^^^^ is given in the inset. 
A clear linear regime can be observed for small V where 
the effective temperatures approach the temperature of 
the leads T^^^^j j j^jT :^ 1. As ^/T increases and the 
non-linear regime sets in, differences between the different 
fixed points are observed. In the Kondo and critical cases, 
all the effective temperatures increase with VjT. In the 
LM fixed point the linear regime where T^^jj^ T ex- 
tends to much larger values of V. This can be intuitively 



understood by the fact that, in this case, the impurity 
spin is mainly interacting with the bosonic bath which is 
insensitive to the chemical potential drop. 

A somehow surprising result is the fact that by replac- 
ing T by Teff,^ the equilibrium universal scaling form of 
the spin-susceptibility is recovered for all cases. Fig. [2]-left 
panel shows that the imaginary part of the susceptibility 
as a function of frequency follows the same equilibrium 
scaling form with T substituted by Teff^^. The static sus- 
ceptibility as a function of Teff,^ (middle panel) also follows 
the same equilibrium universal curve. 

This suggests that Teff,^ is a useful concept to interpret 
the dynamic and static susceptibility out of equilibrium, 
even in the non-linear regime where a single effective tem- 
perature could not be defined. 
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Fig. 2: Left panel - Imaginary part of the spin-susceptibility 
x' {(^) as a function of frequency rescaled by Tgff,;^ com- 
puted for the Kondo = 0), the critical [g — gc) and the 
LM {g = 2.5^c) fixed points for different temperatures and bias 
voltages. The colored points are numerical results obtained for 
different temperatures (see inset caption) and the dashed black 
lines are fits to the equilibrium form. Middle panel - Static 
spin-susceptibility x (^ = 0) as a function of Teff,;^. Points 
with the same color are computed for the same values of T and 
different values oi V 7^ 0. The black lines display the equi- 
librium result. For both quantities x" (^) ^iid x (^ = 0)? the 
equilibrium scaling form holds by replacing T by T^f^^^- Right 
Panel - Comparison between the different effective tempera- 



tures, computed by Eq.(12), as a function of V/T for a fixed 
T. The inset shows the scaling collapse of Teff,;^/T for different 
values of T and V . 

I-V Characteristics. We now turn to the discussion 
of the IV characteristics, see Fig. [3l For V ^ the con- 
ductance increases for low temperatures and approaches 
a constant. Go at zero temperature. Fig. |3] shows the 
behavior of the current for the three regimes for a non- 
equilibrium setup with AT = 0, F ^ 0, computed at 
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different temperatures. Here we define the linear conduc- 
tance per channel as 



G(T) 



^ dJn 

M dV 



For ^ = the current is proportional to the applied voltage 
GqV as long as y,T <C Tk- Outside of the scaling 
regime, i.e. for V^T > Tk^ Jn/y drops rapidly when V or 
T increase. For < g < Jn/V still approaches Go for 
small voltages, however the drop in the conductance arises 
for a smaller energy scale V^T T^{g) where T^{g = 
0)=Tk andT^(^ = ^e) = 0. 

In the LM regime the T = conductance vanishes as the 
impurity effectively decouples form the conduction elec- 
trons. The dependence of the current as T, F ^ 0, in the 
linear and non-linear response regimes develop different 
power law behaviors. The exponents are well fitted by the 
scaling ansatz of Ref. |6 : in the linear response regime 
Jn (X T^-"* V and in the non-linear case Jn cx V^-"*. 

In the critical regime {g = gc) the relation between the 
current and the applied voltage is still linear (j7n = GcV) 
as first found in 6 , however the typical values of the criti- 
cal conductance Gc are much smaller than Gq. The linear 
and nonlinear regimes are characterized by slightly differ- 
ent values of the conductance separated by a crossover. 
In Ref. 6 a similar behavior was reported but a larger 
variation between the linear and non-linear conductance 
was found. This can be explained by the fact that Ref. 6 
considers the current computed using the T-matrix from 
the underlying Anderson model. 

(h) Finite temperature gradient. Fig JI] shows the de- 
pendence of X Te// on the out-of-equilibrium drive 
parametrized by AT/T where T = {Tl + Tr) /2 is the av- 
erage temperature. In this case, for all the critical points, 
the effective temperatures are strongly affected by changes 
in AT/T, since fermionic and bosonic excitations are sus- 
ceptible to a gradient in T across the impurity. The agree- 
ment between the effective temperatures computed for the 
three different quantities is much better than in case (a) 
(where AT = 0,1/ ^ 0). 

As for the setup (a), x as a function of Teff,;^, is found 
to follow the equilibrium scaling form for all regimes. The 
fact that Te// is now strongly renormalized at all fixed 
points shows that the scaling behavior is robust. 

Thermal transport. Due to the particle-hole symmetry, 
AT across the leads does not induce a net particle current 
but there is an energy flow from the hot to the cold lead. 
We consider the fermionic contribution to this energy flow. 



defined in Eq.(ll) As the bosonic contribution to the en- 



ergy current turns out to be vanishingly small, we focus 
on the fermionic part of the energy flow. 

Fig. |5]shows the scaling form of the energy current Je = 
Je,c,L as a function of AT/T for the different regimes. In 
the Kondo regime, for g = 0, J7e/T scales linearly with 
AT/f as long as T < T^. For T > T^, Je/f drops to 
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Fig. 3: Left Panel - Scaling of the particle current Jn as a func- 
tion of V/T for the different fixed points computed for several 
temperatures (the color code follows the one ofFig[2]). Right 
Panel - Jn/V normalized to the conductance of the unitary 
limit Go as a function of temperature (upper panel) and bias 
voltage (lower panel) for different values of g. 



zero. In the hmit AT ^ and T < T^, 



K 



1 1 dJe 



MT dAT 



AT=0 



converges to Kq for g < gc. In the critical regime, K 
remains finite in the AT limit, with Kc < Kq. For the 
LM fixed point K vanishes for f ^ as Je cx ATT 

Conclusion. — In this work we analyzed the concept 
of effective temperatures near local quantum criticality 
focusing on two non-equilibrium setups ((a):]/ ^ and 
AT = and {h):V = and AT ^ 0) and studied the 
nonlinear energy and charge currents in the system. The 
model considered here, the generalized Bose- Fermi Kondo 
model in a dynamical large N limit, can be solved ex- 
actly. We find that for all scaling regimes everywhere in 
the phase diagram and for all considered non-equilibrium 
setups, the equilibrium scaling form of the static and most 
remarkably of the dynamic spin susceptibility can be re- 
covered by utilizing the effective temperature in the equi- 
librium scaling relations, rather then the temperature of 
the leads. The effective temperature as a function of the 
non-equilibrium drive, has been found to be qualitatively 
different in setups (a) and (b). As local observables can be 
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Fig. 4: Effective temperatures for a steady state obtained by 
imposing AT / computed in the asymptotic low frequency 
limit from the FDR of the pseudo-particle propagators and 
susceptibility. 



computed by Wick's theorem, their effective temperature 
is completely determined by the effective temperatures of 
the pseudoparticles. Our results therefore suggest that 
steady-state response near (local) quantum criticality ap- 
pears thermal albeit an effective temperature. The model 
considered here is one example of unconventional quan- 
tum criticality 
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which has recently been discussed in 
the context of holographic duals [26] . Very recently, it has 
been suggested that for quantum critical systems possess- 
ing a gravity dual, out-of-equilibrium current noise can 
appear thermal 
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It would be very interesting to extend the analysis pre- 
sented here to other models of steady state quantum crit- 
icality and explore the extent to wich our findings are 
generic for quantum critical systems. 
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